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METHOD FOR OPTIMIZING THE TRANSMIT SIGNAL IN MULTIPLE 

ANTENNA WIRELESS LINKS 

BACKGROUND OF THE INVENTION 

[0001] The use of multiple transmit and receive antennas can 
enable very large increases in the capacity — defined as the throughput 
or bit rate that can be reliably conveyed — of wireless communication 
systems. For a given set of transmit and receive antennas, the highest 
capacity is attained when the transmit signal in the multiple antenna 
wireless links is optimized. This is straightforward when both the 
transmitter and the receiver can track the instantaneous state of the 
radio channel between them. Such tracking is fairly simple at the 
receiver, but very challenging at the transmitter, particularly in mobile 
systems, because of the fast-changing nature of the channel. Namely, 
by the time the receiver can report the instantaneous state of the 
channel to the transmitter such that the transmitter can adapt to the 
instantaneous state, the instantaneous state has changed. Hence, it is 
often the case that the transmitter must operate without up-to-date 
reliable information about the instantaneous response of the channel, 
and sub-optimal covariance and thus capacity is achieved. 
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SUMMARY OF THE INVENTION 
[0002] The present invention provides a method of optimizing 
transmit signal for a multiple antenna transmitter using a statistical 
distribution of the channel. As stated above, the receiver has 
knowledge of the instantaneous state of the channel. This is typically 
because the transmitter transmits pilot signals known a priori to the 
receiver, which may then derive the instantaneous state. As further 
stated above, the transmitter lacks knowledge of this instantaneous 
state. Statistical information about the channel, however, is available 
because the channel is short-term stationary and thus its statistics 
remain unchanged for an amount of time that is hundreds or even 
thousands of times longer than the amount of time over which its 
instantaneous state varies.. The distribution is the set of possible 
instantaneous states and the probability associated with each one , 
where each probability indicates the probability that the associated 
instantaneous state currently characterizes the channel. 
[0003] In one exemplary embodiment, spatial steering of parallel 
symbol streams and a transmit power of each symbol stream is 
determined based on a statistical distribution of a channel. For 
example, respective steering vectors for the parallel symbol streams 
aire determined as the eigenvectors of the expected distribution of the 
channel matrix multiplied by its transpose conjugate. 
[0004] In an exemplary embodiment, the transmission power for 
each symbol stream is determined based on a previously determined 
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transmission power for the symbol stream and the statistical 
distribution of the channel. In this manner, the determined 
transmission power for each symbol stream finds and tracks the 
optimal transmission power for each symbol stream. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0005] The present invention will become more fully understood 
from the detailed description given herein below and the 
accompanying drawings, wherein like elements are represented by like 
reference numerals, which are given by way of illustration only and 
thus are not limiting of the present invention and wherein: 
[0006] Fig. 1 illustrates an exemplary structure of a transmitter to 
which the method of the present invention may be applied. 

DETAILED DESCRIPTION OF THE EXEMPLARY EMBODIMENTS 
[0007] A multiple antenna wireless system having M transmit and N 
receive antennas may be modeled at baseband as: 

y = Hx + n 

[0008] where x is an M-dimensional vector representing the 
transmit signal that we want to optimize, H is an NxM matrix 
representing the channel between each transmit and receive antenna, 
y is an N- dimensional vector representing the received signal and n is 
another N-dimensional vector that represents the noise. The noise is 
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modeled as Gaussian and white, but the generalization to non- white 
noises is straightforward. The channel, in turn, is considered 
frequency-flat over the bandwidth of x but, again, the extension to 
frequency-selective channels is straightforward by one of several well- 
known methods (such as breaking x into a number of disjoint narrow 
sub-signals, equalizing y to flatten the channel, etc). With a 
frequency-flat channel, the entries of x, y, H and n are scalars, and, in 
general, are complex to accommodate in-phase and quadrature 
components. 

[0009] As stated above, the receiver has knowledge of the matrix H 
at all times. This is typically achieved by the transmitter transmitting 
pilot signals known a priori to the receiver, from which H can be 
derived. As further stated above, the transmitter lacks knowledge of 
the matrix H. Statistical information about H, however, is available 
because H is short-term stationary and thus its statistics remain 
unchanged for an amount of time that is hundreds or even thousands 
of times longer than the amount of time over which H itself varies. The 
instantaneous response of the channel can be used to derive the 
distribution of the matrix H in any well-known statistical manner. The 
distribution of H is the set of possible values H can take and the 
probability associated with each one, where each probability indicates 
the probability that the associated instantaneous value corresponds to 
the actual H. 
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[OOIO] This distribution of H changes slowly over time and thus the 
transmitter must also slowly vary its settings accordingly. We shall 
indicate by T the duration of time over which the covariance of the 
generated signal remains unchanged and we index by a counter k the 
succession of such intervals. 

[0011] Also, we denote by P the power available for transmission 
and thus we have that 

E[x H x] = E[Tr{xx H }] = P 

where () H indicates the transpose conjugate of a vector or matrix, Trf} 
indicates the trace of a matrix, and E[ ] indicates expectation, i.e. 
average over the distribution of some quantity. Finally, a 2 indicates 
the arithmetic mean of the noise variances at the N receive antennas. 

[0012] Fig. 1 illustrates an exemplary structure of a transmitter to 
which the method of the present invention may be applied. This is a 
well-known structure that enables the generation of a signal with 
arbitrary structure, and will provide a context for better understanding 
the method of the present invention. As shown, a space-time encoders 
and modulators block receives a stream of bits (referred to as source 
bits) that are to be conveyed to a receiver reliably. These bits are 
encoded and then mapped onto M complex modulated symbol streams 
in a well-known manner. The amount of information embedded onto 
each symbol in a symbol stream is determined by a rate R m (k), 
m=l,...,M provided to the space-time encoders and modulators block 
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10 for each respective symbol stream. At each interval k there are 
thus M such rates. Some of these rates may be zero at times 
indicating that the corresponding symbol will not be used during that 
interval and hence no information should be embedded therein. These 
rates may be supplied to the transmitter by the receiver or they may 
be computed directly at the transmitter (in which case the structure of 
the receiver must be known therein). Because the operation and 
structure of the space- time encoders and modulators block 10 is well- 
known in the art, this block will not be described in detail for the sake 
of brevity. 

[0013] Each of the M symbol streams (also interchangeably referred 
to herein simply as "symbols") outputted by the space-time encoders 
and modulators block 10 is weighted by an associated distinct 

amplitude ^p m {k) , m=l M. As shown in Fig. 1, each of the M 

symbol streams is multiplied by ^jp m (k) , m=l,...,Mat multipliers 12. 
At every interval fc, these amplitudes must satisfy 

so that the power budget P is not exceeded. Namely, the sum of the 
transmission powers for the symbols streams can not exceed the 
power P available at the transmitter. Thus, p m (k)/P indicates the 
fraction of the total transmitted power P that is allocated to the m-th 
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symbol during interval k. For convenience, we can assemble the 
transmission powers onto the following diagonal matrix: 

>,(*) 0 ••• 0 
P0d= ° Pz(k) Pl{k) 1 

; o -. "o 

.0 ••• 0 p M (k) 

[0014] After being allocated its power, the m-th symbol is multiplied 
by a vector [Vi.m(fc) V2. m (k) ... V M ,m{k)] which steers it spatially. As 
shown in Fig. 1, a set of M multipliers 14 is associated with each 
symbol stream, and multiplies the symbol stream by the associated 
steering vector. These M steering vectors, each of which has M 
complex coefficients, may be assembled into an MxM matrix as shown 
below. 



V(k) = 



2,2 (*) 



V u (k) V h2 (k) 

v w (*) v. 

': 0 



V M . M (ft) 



[0015] The matrix V, in an exemplary embodiment, is unitary, i.e. 
such that W H =I, where I denotes the identity matrix. Consequently: 
a) Each of the steering vectors (each column of V) has unit norm, 
that is, 



M 



£K.m| =1 f° r ™ = \,...,M 
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and therefore the power allocation P is not altered. Each symbol 
retains its power through the steering process. 

b) The M steering vectors (columns of V) are mutually orthogonal 
and thus they do not interfere with each other. 

[0016] The last step in the steering process is the combining. The 
sum of the j-th entry of each of the M steered symbols results in the 
signal that will be transmitted by the j-th antenna, for j=l,...,M. As 
shown in Fig. 1, an adder 16 associated with each of the M antennas 
receives the associated entry in each of the M steered symbols, and 
combines them to form a signal for transmission. 
[0017] Before being radiated, the M signals composed by the 
steering process are filtered and upconverted onto their appropriate 
frequency by an associated RF block 18. Finally, the resulting RF 
signals are fed onto the M transmit antennas 20. 
[0018] Next, an exemplary embodiment of detemiining the 
transmission powers and steering vectors to optimize the transmit 
signal in order to maximize the capacity of the wireless links, 
according to the present invention will be described. As discussed 
above, the transmission powers may be characterized as a matrix P 
(whose diagonal entries describe the power allocation) and the steering 
vectors may be characterized as a matrix V (whose columns are the 
steering vectors). The covariance of the transmitted signal equals 
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VFV H , the application of the invention results in a signal whose 
covariance is optimum in a capacity-sense. 

[0019] At interval k, the columns of V(k) are computed to be the 
eigenvectors of ElH^H], where the expectation is computed over the 
distribution of H during the interval. 

[0020] The computation of P is performed via a recursive algorithm 
which updates it at every interval based on (0 the distribution of H 
during the interval, and (u) the value of P at the previous interval. In 
acquisition mode, P is initialized arbitrarily and then the optimum 
powers are found after a few iterations. In steady state, the algorithm 
tracks the optimum set of powers as the channel distribution slowly 
varies. 

[0021] In order to formulate the algorithm, it is useful to 
introduce — as an auxiliary quantity — the mean-square error that 
would be experienced when estimating with a linear minimum mean- 
square error filter the m-th transmitted signal (i.e., the signal steered 
by the m-th column of V). Defining 



where H = HV while H_ m and P_ m O) indicate the corresponding 
matrices without the m-th column, such mean-square error is given by 
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MSE(k) = 



1 + ^H"B H 



with H m denoting the m-th column of H . 
The algorithm is as follows. At interval Jc+1: 
1) The first step is to compute, for m=l,...,M, 



P m (£ + l) = max 



E[Tr[B^}hpi^^-N 



with the expectations taken with respect to the distribution of H 
during interval k 

2) The second and final step yields the new value of P as 



P(k+1) = 



M 



m=l 



P\ 0 0 
o P 2 o 
0 0 '•. 



0 
0 
0 



o 0 0 p' M 



[0022] Namely, the second step normalizes the transmission powers 
determined in the first step such the available power at the 
transmitter is not exceeded. 

[0023] In one exemplary embodiment of the present invention, the 
method described above is implemented as software in the processor 
of the receiver. The receiver, at each interval k, would convey the 
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resulting V(Jc) and P(fc) to the transmitter. In an alternative 
embodiment, the method is implemented as software at the 
transmitter with the receiver sending the transmitter the necessary 
information to determine the matrices V(Jc) and P(fc). As further 
alternatives, the method according to the present invention may be 
implemented in hardwired circuits, firmware, FPGAs, etc. 
[0024] The invention being thus described, it will be obvious that 
the same may be varied in many ways. Such variations are not to be 
regarded as a departure from the spirit and scope of the invention, 
and all such modifications as would be obvious to one skilled in the 
art are intended to be included within the scope of the following 
claims. 
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